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ABSTRACT 

Transport of angular momentum is one of the thrust areas of astrophysical flows. 
Instabilities and hence the turbulence generated by it has been invoked to understand 
its role in angular momentum transport in hydrodynamic and magnetohydrodynamic 
regime. Investigation of an unexplored region described by the parameter space fl^ < 
0, Oe being the epicyclic frequency, resulted in a powerful magnetohydrodynamic 
instability. The growth rate of this instability is rather large and is found to depend 
on the wavenumber. 
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1 INTRODUCTION 

Transport of angular momentum i n various as tr ophys - 
ical flows is of great im portance (|Lvnden-Belll ([l96^; 
IShakura fc SvunyaevI l| 19731 )) in accretion theory. Nature of 
accretion process is largely understood through theoretical 
and computational studies. In astrophysical flows, outward 
transport of angular momentum of accreting material is the 
basic requirement to sustain the accretion disk. This angu- 
lar momentum transport can be sustained by some turbulent 
mechanism operating in accretion disk. The origin of turbu- 
lence is unclear. Turbu lent viscosity (g-visco s ity) m odel by 
Shakura and Sunayev (|Shakura fc SvunyaevI (| 19731 )) IS one 
of such mechanisms. Differential rotation and the presence 
of magnetic field plays an important role in the dynamics of 
instability and hence in efficient angular momentum trans- 

f ort. In early 1990s, Balbus & Hawley (jBalbus fc HawlevI 
199ll )) found that the subthermal fields in the presence of 
differential rotation can cause Magneto-Rotational Instabil- 
ity (MRI) , which can provide the physical basis for the trans- 
port of angular momentum in accretion disk. Although the 
MRI requires only a small amount of ionization to work, 
it cease to exist when the strength of magnetic field in- 
crease s. Paper by Bonanno and Urpin l|Bonanno fc UrpinI 
l|2006t )) (henceforth paper I) is an attempt to understand 
the transport of angular momentum by studying the insta- 
bility arising in differentially rotating compressible flows in 
the presence of high magnetic field. The magnetic field had 
nonvanishing radial and azimuthal components. They found 
a new instability which survived for any value of magnetic 
field unlike MRI which survives only in the weak field limit. 
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The maximum growth rate of the reported instability was 
~ where Q is the rotation frequency. 

This work investigates the magnetohydrodynamic in- 
stability in an relatively unexplored region in accretion disk 
described by < 0, Qe being the epicyclic frequency. 

For brevity, we explain the physical meaning of epicyclic 
frequency. Square of Epicyclic frequency is defined as Q,^ — 
4f2^(l -I- 2n^)i s being the radial distance from the rota- 
tion axis, Q being the angular velocity of the fluid element 
in the flow. The physical meaning of the epicyclic frequency 
is explained as follows: Orbits of stars in a galaxy are usu- 
ally not exact circles and always have a certain amount of 
randomness. If one of the perturbed orbits is observed on a 
framework rotating with the mean angular velocity around 
the cen t er of the galaxy, it is known to be a small circle 
iBertinl l| 19961 )). The frequency of such small cyclic motions 
that originate from randomness is called the epicyclic fre- 
quency. 

Negative epicyclic frequency can be understood as fol- 
lows: Square of epicyclic frequency is defined as 



(1) 



specific angular momentum. An inward increase in the spe- 
cific angular momentum (i.e, outward decrease in angular 
momentum) implies that the square of epicyclic frequency 
is negative. A fluid element displaced outward under conser- 
vation of angular momentum has a larger angular momen- 
tum compared with that of the surrounding media. In this 
case, the centrifugal force acting on the displaced element is 
stronger than that acting on the surroundi ng me d ia, an d the 
element goes further outward (unstable) (|Katd (|2007l )). In 
literature, it is known as a Reyliegh criterion for instablity. 
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It is well known that in a keplerian disk, the specific 
angular momentum increases with radius, i.e., > 0. Ac- 
cording to the Rayleigh criterion, such a flow is stable. Gen- 
erally in accretion disk, the region with £7^ > do not exist 
but su c h regions are realized loca l ly in the models bvlHonmal 
j 19961 ) : lManmoto fc Katol (|2000|): iKato fc Manmot ol l|2000D . 
A recent work by Gracia et al.l 1 20031 ) has numerically shown 
that the transition region between ADAF and outer stan- 
dard disk is time-dependent having Q,^ < 0. They modelled 
the time-evolution of bimodal accretion discs where an outer 
cooling-dominated disc matches an inner ADAF type flow 
through a narrow transition region. The behaviour of the 
instability in such regions is explored. 

This paper is orgnized as follows: The section 2 investi- 
gates the the growth rate of instability in the new region of 
parameter space. The section 3 deals with the Results and 
Discussions. Finally, we conclude by summarizing the new 
findings. 



2 THE INSTABILITY CRITERIA 

Paper I considered an axisymmetric different ailly rotating 
system in the presence of a magnetic field. A cylinderical 
coordinate system (s, ip, z) with s being the radial distance 
from the rotation axis was constructed. Unperturbed system 
was decscibed by (fr, Wz) = (0, sr2,0). Furthermore, SI 
where being the angular velocity of the astrophysical flow, 
was taken to be approximately a function of s alone; i.e., 
Q, = Q.{s). The isothermal flow for a compressible fluid was 
described by the following MHD equations 



(iT- V)« = -f-ff-l- --!-(V X B) X B, (2) 

P 4:TVP 



p + \7-{pv)=0, (3) 
p + V ■ Vp + 7pV ■ V — 0, (4) 

V X (iTx B)-hJ?V X (V X B) =0, (5) 

V ■ B = 0. (6) 

Where p and v are the density and fluid velocity, respec- 
tively; p is the gas pressure; g is gravity; B is the magnetic 
field, rj is the magnetic diffusivity, and 7 is the adiabatic in- 
dex. Magnetic field has non vanishing radial and azimuthal 
components. 

Quasistationary state for a differentially rotating fiow 
with magnetic field was considered. Axisymmetric Eulerian 
perturbations with space time dependence oc exp{at — ik-r) 
where k = (fcs,0, fcz) were introduced in the unperturbed 
accretion disk and the dispersion relation was obtained for 
a special case of k.B — after neglecting the Ohmic dissi- 
pation in the induction equation. 

We take the dispersion relation of paper I for determin- 
ing the stability of axisymmetric short wavelength perturba- 
tions . The dispersion relation (Eq. (16) of paper I) is given 
by the 

+ cr^(cJo + f^e) + f^^t^sn + cr^figLJo + p-^^e^^Bn = (7) 
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Equation ([YJI is a polynomial ol degree tive, so 
five non trivial roots exist. We appl y Routh-Hurwitz 
method (|Aleksandrov fc LaurentievI (|l985h) (see appendix) 
to find the regions describing instabilities. This methos 
has been applied in the fi e ld of astrop h ysics ve r y fre- 
queri tl y jBonanno fc UrpinI ||2006|. |2007|): balbusl (|2000l . 
l200j ): lUrpin fc Riidigerl l|2005h : lUrpinI (|2003l )) Instabilitv 
exist if any of the conditions written below is satisfied. 



pilluj%Q < , c^ln > , {u}%nf > 0. 



(8) 



From the above inequalities, it is interesting to note that 
the instability will exist for fi^ < also. We will investigate 
this region for instability. 

Polynomial given by Eq. in < region can be 
written as 



Let us explore the behavour of dispersion relation for 
incompressible flow. Since incompressibility corresponds to 
00, Equation ^ becomes 
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(10) 



It is to be noted that only those perturbations are con- 
sidered here in which the wavevector is perpendicular to the 
magnetic field, i.e. fc • _B = 0. It is interesting to note that in 
the compressible limits, the magnetic field, irrespective of its 
strength, do not play any role. The hydromagnetic instabil- 
ity cease to exist and the fluid behaves as a Rayleigh unsta- 
ble hydrodynamic entity. This indicates that in the presence 
of magnetic field, the compressibility of fluid plays a major 
role in deciding the occurence of hydromagnetic instability 
for the above dicussed magenetic field configuration. This 
behaviour, i.e., the dependence of hydromagnetic instability 
on the nature of fluid (compressible/ incompressible) can be 
understood by recalling that the magnetorotational instabil- 
ity is a property of incompressible fluid and it cease to exist 
for the compressible flow. 

When we consider the non magnetic case, the dispersion 
relation yields Reyliegh unstable modes. 

Let us investigate the growth rate of the hydromagnetic 
instability. To calculate the growth rate of this instability in 
this unexplored region, < 0, it is convenient to introduce 
the dimesionless quantities 



r = 
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m 



c = 
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X — ks 



(11) 



The polynomial given by Eq. Q in the case of < can 
be written as 



1) +r\a;^ 



r^^x^ - ^l(x^ = (12) 



This equation is solved numerically (see iPress et al.l 

l|l992l ) for details) by computing the eigen values of the ma- 
trix whose charesteristic polynomial is given by equation ((9| 
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Figure 1. The dependence of real and imaginary parts of T on Figure 3. The dependence of real parts of T on for fj, = 0.3, 

for ^J, = 0.3, e, = 0.1, C = 0.1 f = 0.1 and C = 0.5, 1.0 



30 
20 
10 

¥ ° 
--10 

-20 

-30 

-40 



100 200 300 400 5^0 600 700 800 900 1000 



Figure 2. The dependence of real and imaginary parts of F on 
for n = 0.3, ( = 0.5, C = 0.1 
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Figure 4. The dependence of real F on x^ for fj, ■ 
and c; = -0.1 



0.3, C = 0.1 



for different values of parameter ^, ^ and We get three 
real roots and one pair of complex root unlike paper I where 
two pairs of complex conjugate roots ware obtained. 

Fig. 1, shows the dependence of real and imaginary 
parts of r on X for jj, = 0.3, ^ = 0.1, ^ = 0.1. Two real 
roots are obtained when ^ — 0.3, £, — 0.1 and ( — 0.1 whose 
growth rates are ~ 0.5Qe but when value of ^ is increased 
(Fig. 2), growth rate of one of the roots with positive real 
part deminishes and for any further increase in the value of 
^, this root decays very fast and approaches zero but the 
growth rate of the other root remains around ~ O.Sile. For 
^ >> ^, growth rate of one of the roots approaches zero as 
in paper I but the growth rate of the other root with posi- 
tive real part remains F ~ 0.5 and is almost independent of 
the wavenumber beyond a certain value. It appears that the 
root in paper I which diminishes with increasing value of ^ 
is similar to the root which diminishes in our calculations 
also (it is clear if we compare Fig, 1 of our paper with the 
Fig. 1 of paper I) but the second surviving positive root in 
our calculation is possibly a contribution oi Q,^ < 0, which 
was neglected in the paper L This is qualitatively clear be- 
cause ^ » corresponds to incompressible limit where the 
reported instability of paper I diminishes and only contri- 
bution from f2g < remains. 

But a very interesting result is obtained when we cal- 
culate the growth rate of instability for ( > Value of 
^ was fixed to 0.1 and the parameter ^ was varied with 
C = (0.5,1.0,1.5,2.0,3.0). It is observed that the root in 
this case is growing faster than the previous cases and this 
growth is not independent of the wavenumber. Fig. 3, clearly 



shows this pattern. Interestingly the growth rate of the in- 
stability is as high as 5.9n^ for =900 at C = 3.0. 

Undoubtedly this is a very strong instability and much 
more powerful than what was observed in the paper I. The 
maximum growth rate reported in Paper I was ~ Qe for 
J5 200. The instability in our case is very important on 
two counts, i) Growth rate is very high (~ 5.9f2e) for x^ — 
900 at = 3) ii) unlike the previous report, it depends on 
the wavelength of perturbation. 

In Fig. 4, we plot the dependence of Re (F) on x for 
H = 0.3, ^ = 0.1 and ^ = —0.1. It is very clear that the 
change of sign of reverses the behaviour of roots (see Fig. 
1, and Fig. 4.). We do not plot the growth of imaginary part 
of root as we are interested in the instability. 



3 RESULTS AND DISCUSSION 

We revisit the problem of hydromagnetic instability in dif- 
ferentially rotating compressible flows and analyse the insta- 
bility. The behaviour of hydromagnetic instability is investi- 
gated in the region described by Q,^ < 0. Though 0^ < is 
generally not satisfied in accretion dis k but there ar e poss i- 
bilities that < is lo c allv satisfied |Gracia et al.l |2003): 
iManmoto fc Katd l|2000l ): lKato fc Manmotd l|2000l ): lH'onmal 
in the accretion disks. The present work will be ap- 
plicable in those regions (Kato (2Q0S))- 

We discover a range of parameter space described by ^, 
^ and fj, where the instability grows very fast, much faster 
than the reported growth rate of the instability in paper I. 
We have shown the growth rate of the instability in Fig. 3, 
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for two representative values of ^ for fixed Growth rate 
is ~ l.Sfle, ~ 2.0r2e and ~ 2.fr2e at x'^ = 300, 600 and 900 
respectively for ^ = 0.5 and the growth rate is ~ 2.6il,e and 
~ 3.0f2e and ~ 3.4f2e respectively for same values of x'^ for 
= 1.0. The growth rate increases further with increasing 
value of for fixed values of ^. It is surprising to observe that 
the instability grows very fast and the growth rate, unlike 
in paper I is not independent of wavelength of perturbation, 
Infact, the instability grows with the increasing value of x^. 

In the case of ^ >> i^, which corresponds to the incom- 
pressible limit, one of the roots diminishes as was in paper 
I but the second root with postive real part survives with 
growth rate 0.547732f2e. Since incompressibility corresponds 
to tt;o ^ oOi Equation ((9)1 becomes 

(J (cr^ - Atl^^el) = (13) 

It is imperative to see that in the case of incompressible 
flow, Eq.® reduces to hydrodynamic relation. Since jj, = 
0.3, growth rate comes to ,/fIQ,e, i.e., .0547732. Hence, it 
is clear that the root which survives in the incompressible 
limit corresponds to the root satisfying Rayleigh criterion of 
instability in hydrodynamics. 



4 CONSLUSION 

This paper considers a generalization of the well known 
Rayleigh instability for the case of compressible magnetic 
flows. The growth rate of the instability in the new region 
f2g < is rather large. Moreover, the instability in our case 
also depends on the wave number. High growth rate and de- 
pendence on wavenumber makes this instability very inter- 
esting. Since the instability exists even in the high magnetic 
field where MRI is supressed, it will be interesting to study 
the effect of this instability in angular momentum transport. 
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APPENDIX A: HURWITZ METHOD 

Let us consider an 5*'' order polynomial P(x) 

P{x) — a^x^ + a^x"^ + a-iX"^ + a2X^ + + ao (Al) 

Then Hurwitz theorem states that the above polynomial will 
be unstable if any of th e following inequalities is satified 
l|Urpin fc Riidigej l|2005^ ). 

ao < , 

Ai = a4a3 — a2 < , 

A2 = 02(0403 — 02) — 04(0401 — oo) < , 

A3 = (0401 — 00) [02(0403 — 02) — 04(0401 — Oo)] — 

-00(0403 - 02)^ < , (A2) 



